Nonaxial-octupole effect in superheavy nuclei 
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The triaxial-octupole Y32 correlation in atomic nuclei has long been expected to exist but experimental ev- 
idence has not been clear. We find, in order to explain the very low-lying 2~ bands in the transfermium mass 
region, that this exotic effect may manifest itself in superheavy elements. Favorable conditions for producing 
triaxial-octupole correlations are shown to be present in the deformed single-particle spectrum, which is fur- 
ther supported by quantitative Reflection Asymmetric Shell Model calculations. It is predicted that the strong 
nonaxial-octupole effect may persist up to the element 108. Our result thus represents the first concrete example 
of spontaneous breaking of both axial and reflection symmetries in the heaviest nuclear systems. 

PACS numbers: 21.10.Re, 21.60.Cs, 27.90.+b 



The issue that underlies the present investigation is whether 
some heaviest atomic nuclei can be characterized as having 
an intrinsic mean field with nonaxial-octupole component, in 
connection with the tetrahedral symmetry which is quite com- 
mon in molecular, metallic clusters, and some other quantum 
objects. The tetrahedral symmetry is a direct consequence 
of the point group and corresponds to the invariance under 
transformation of the group T^ which has two one- and one 
four-dimensional irreducible representation. In molecular and 
metallic clusters, the tetrahedral symmetry is determined by 
mutual geometric arrangement of the constituent ions. In 
atomic nuclei, which are finite many-body systems with the 
strong interaction, the situation is more complex. However, it 
is usually understood that nuclear shapes are governed by the 
shell effects. The tetrahedral instability in finite many-fermion 
systems was suggested in 1991 QJJ] . Tetrahedral shape in 
nuclei was predicted with more realistic model calculations 
I2,3,fli. 



For description of nuclear shape, it has been proven to be 
extremely useful to parameterize the surface R(9,(j)) through 
its expansion in the spherical harmonic function Y^(0,0) , 
where 9 and (j) are the rotation angles [6]. The coefficients 
of the expansion , a^, are related to deformation parameters. 
The quadrupole deformation (A = 2) is the leading deforma- 
tion effect known almost for all the nuclei throughout the nu- 
clear chart. The next one, the octupole deformation (A = 3), is 
related to a spontaneous breaking of reflection symmetry, or 
parity, and is known in some nuclei. The tetrahedral symmetry 
is realized at the first order through the nonaxial-octupole Y32 
deformation J3, Si • Almost pure Y32 shape was suggested by 
the two level model study of finite Fermion systemsfl. Cal- 
culations of possible nonaxial pear-like ground-state deforma- 
tions in the mass region 220 < A < 230 has been reported ifioll . 
Both axially- and nonaxially-enhanced octupole correlations 
have been suggested for the A — 90 mass region [11] and in 
the superdeformed Hg nuclei II 1211 . The Skyrme HFB calcu- 
lation has predicted that some nuclei can be extremely soft 
against the Y32 tetrahedral deformation Il3l[l4ll . 



One may thus conclude that there have been many theo- 
retical discussions about this exotic type of nuclear deforma- 
tion; however, experimental evidence has not been very clear. 
This may imply that the tetrahedral symmetry and the pure 
Y32 shape manifest themselves in nuclear physics in a very 
complex way. The octupole mode, as a smaller deformation 
effect, has to compete with the dominant quadrupole mode, 
and as a result, the appearance of the tetrahedral symmetry in 
all known mass regions is somewhat obscure. The question 
is whether one can find evidences in the heaviest mass region 
which is still less known. It is legitimate to ask this question 
because a strong octupole effect has been well known to ex- 
ist in the actinide region B15I1 which is the closest neighboring 
mass region to the heaviest nuclei. The present study will try 
to answer this question. Some experimental data, though very 
limited and incomplete, have already indicated such a possi- 
bility. These data show that very low-lying negative parity 
bands with the sequence of I 71 = 2~,3~,4~,5~, • • • appear in 
the spectrum of some very heavy nuclei (see discussions be- 
low). 

Exploring the predicted stable island of superheavy ele- 
ments is the current goal in nuclear science. In the past few 
years, progress has been made in synthesis of new elements 



(for review, see Refs. |3 17, 1_8 , 19]). On the other hand 



the heaviest nuclei for which detailed spectroscopic measure- 
ments can be performed are the isotopes with proton num- 
ber Z - 100 and neutron number A? ~ 150 EH 0. The 
significance of studying the deformed nuclei in this transfer- 
mium region is that it can provide an indirect way to access the 
single-particle states of the next spherical shells because some 
of these orbitals are strongly down-sloping, and thus can come 
close to the Fermi surface in the deformed region [23]. 

There has been experimental observation of low-energy 7T 
bands in the mass-250 nuclei. For example, for the N = 150 
isotonic chain, early data showed a 2T band in 246 Cm with an 
bandhead energy E(2~) = 0.842 MeV JH, and in 248 Cf with 
a much lower bandhead E(2~) = 0.592 MeV [27]. The recent 
data have added new evidences for 250 Fm with E(2~) = 0.879 
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FIG. 1: Experimental energies of the 1 _ (solid symbols) and 2~ 
(open symbols) bandheads of the isotopes Ra, Th, U, Pu, Cm, Cf, 
Fm, and No. Data for Ra, Th, U, and Pu are taken from Refs. 
(HE!, and for Cm, Cf, Fm, and No from Refs. (U, HI, H|, HI- 
The solid line and dash line are the polynomial fits for the 1~ and 
2~ bandhead energies, respectively. Note that the two curves cross 
at A = 242. 

MeV J3 and for 252 No with E(2~) = 0.929 MeV 0]. We 
find it difficult to interpret these bands as 2-quasiparticle (qp) 
excitations in a reflection-symmetric well because of their low 
bandhead energy, particularly for 248 Cf. The reason is sim- 
ple. To be excited, a 2-qp state must have an energy greater 
than two times of the BCS pairing gap, and the pairing gaps 
in this mass region are typically ~ 0.5 MeV for neutrons and 
~ 0.8 MeV for protons, which are the needed amount to cor- 
rectly reproduce the rotational behavior, e.g. the moment of 
inertia. Therefore, any (relatively) pure 2-qp state should lie 
higher than 1 MeV above the ground state. For the observed 
low-energy 2~ bands, there must be strong correlations that 
sensitively act on the negative-parity configurations to push 
them down. The primary candidate of the correlations is of 
the octupole type. 

The purpose of the present article is to investigate the role of 
the nonaxial-octupole correlations in collective excitations of 
transfermium and superheavy nuclei. We first analyze the sys- 
tematical trend of experimental low-lying spectroscopy data, 
associated with axial-octupole and nonaxial-octupole modes 
and ranging from the A = 220 (actinide) to the 250 (trans- 
fermium) mass region. This analysis convincingly shows the 
importance of introducing nonaxial-octupole correlations in 
the description of low-lying 2~ bands in the A = 250 mass re- 
gion. Second, we generalize the Reflection Asymmetric Shell 
Model (RASM) 1 33] to include the triaxial-octupole degree 
of freedom, which is applied to describe quantitatively the 2~ 
bands. Based on these results, we interpret these bands in the 
heaviest nuclear systems as having connection to the predicted 
tetrahedral symmetry. 

Experimental energies of the lowest l - and lowest 2~ 
bandheads of Ra, Th, U, Pu, Cm, Cf, Fm, and No are plotted 
as functions of mass number in Fig. 1. To qualify a "band", 
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FIG. 2: (Color online) Proton Nilsson diagram showing single- 
particle orbitals as deformation £2 varies (left part), and as functions 
of £32 with fixed £2 = 0.24 (right part). In the left part, orbitals with 
positive (negative) parity are shown by solid (dashed) curves. 



we use a criterion which requires that for each band, at least 
three lowest states be assigned experimentally. As a general 
tendency, it is seen from Fig. 1 that the 1 ~~ bandhead energies 
increase with increasing mass number, while the 2~ ones be- 
have like a gradually decreasing function of mass number. In 
the Ra-Pu region, the 1 ~ bands lie very low in energy while 
the 2~ bands are quite high. This reflects the well-known fact 
that in Ra-Pu nuclei, the axial-octupole correlation, associ- 
ated with the Y30 shape, is important. However, the situation 
becomes opposite in the Cm-Cf-Fm-No region where the 2~ 
bandhead energies drop to lower values. The occurrence of 
low-lying 2~ bands indicates that the nonaxial-octupole effect 
becomes important in the transfermium region. 

The polynomial fits for the l - (solid curve) and 2~ (dash 
curve) bandhead energies across each other at A = 242, as 
seen in Fig. 1. This suggests a spontaneous breaking of 
the axial-octupole symmetry around this mass number. More 
precisely, the dominating component of octupole correlation 
changes from the axial Y30 at A < 242 to the nonaxial Y32 
at A > 242. In studies of the nuclear reflection asymmetry in 
actinide nuclei with A < 242, one usually neglects, as a good 
approximation, the nonaxial-octupole deformation and adopts 
only the axial Y30 component in many model calculations for 
the low-lying alternative parity bands (see, for example, Ref. 
l3oTl ). Following the same philosophy, one should apply the 
nonaxial-octupole Y32 deformed mean field as a useful start- 
ing point to model the low-lying 2~ bands in transfermium 
and superheavy nuclei, where the axial-octupole Y30 defor- 
mation component can be neglected as a good approximation. 
Of course, the quadrupole deformation that is superposed with 
octupole deformation must be included to describe the exper- 
imental moments of inertia in both actinide and transfermium 
nuclei. 

In an effective treatment with separable forces, the octupole 
effect is described by the long-range octupole-octupole (O- 
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O) interactions between nucleons, which are produced by the 
matrix elements of Yt,^ between single-particle states with 
Ay = A/ = 3 (j and I are the total angular momentum and or- 
bit angular momentum of single-particles, respectively). The 
condition is found satisfied in several mass regions in the 
nuclear chart, and particularly in those transitional regions 
between spherical and well-deformed shapes. Thus, if the 
Fermi surface of a nucleus lies close to a pair of orbitals with 
Aj = Al = 3, the maximum octupole effect is expected. To see 
whether octupole correlations are present in the heavy mass 
region, we plot in Fig. 2 the deformed single-particle states 
calculated with the Nilsson potential. We note that this proton 
single-particle diagram created by using the Nilsson parame- 
ters of Bengtsson and Ragnarsson [31] is very similar to the 
diagram calculated with the Wood-Saxon potential 112311 . The 
latter diagram has recently been used to assign experimental 
band structure in the odd-proton nucleus 249 Bk 03211 . In the left 
part of Fig. 2, proton Nilsson diagram is plotted as a function 
of the quadrupole deformation £2. The calculation predicts a 
sizable spherical energy gap at proton number 1 14. It is par- 
ticularly notable that below the gap, the positive-parity state 
'13/2 (with / = 6 and j = 13/2) and the negative-parity one 
fl /2 (with / = 3 and j — 7 /2) appear to be nearly degenerate 
at the spherical limit. These two nearly-degenerate states fulfil 
precisely the A j = Al = 3 condition for the octupole correla- 
tion. We note that the near-degeneracy of the proton 113/2 and 
J7/2 states is also a result of the Wood-Saxon single-particle 
diagram 112.311 . 

A careful inspection on the left part of Fig. 2 leads to an- 
other intersting discovery: When going away from the spher- 
ical limit and moving to the deformed region, one can clearly 
see the pairs of orbitals with opposite parities that are very 
close in energy and stay nearly parallel as £2 varies. These are 
the pairs {[633]7/2, [521]3/2} and {[624]9/2, [512]5/2}. The 
two orbitals in each pair satisfy the Aj = Al = 3 condition, 
and in addition, are governed by AK = 2 (K is the projection 
of j on the symmetry axis of a deformed nucleus). Thus these 
orbitals are connected by large matrix elements of F32, imply- 
ing a strong triaxial-octupole correlation. The nuclei, whose 
proton Fermi surfaces fall in the proximity of these orbitals, 
may exhibit an £32 shape. The inserted figure at the upper- 
right corner of Fig. 2 illustrates the most favorable shape of 
a nucleus with triaxial-octupole deformation superposed with 
quadrupole deformation. Due to the strong octupole correla- 
tion the collective negative-parity states (in this case the se- 
quence with the bandhead I % = 2~) can be pushed down sig- 
nificantly. 

It may possibly be true that the strong T32 correlation, with 
the occurrence of very low-lying 2~ bands, is a special prop- 
erty in superheavy nuclei. The strict conditions for this to 
happen are the presence of nearly-degenerate pair(s) of single- 
particle states at £2 = that differ in the quantum numbers / 
and j by 3, and due to this degeneracy at the spherical limit, 
the presence of closely-lying, opposite-parity, sub-orbitals for 
£2 > that differ in the quantum number K by 2. 

In the right part of Fig. 2, the diagram is obtained by includ- 



ing an additional 132 term in the Nilsson Hamiltonian vary- 
ing with the octupole deformation parameter £32, while the 
quadrupole deformation is kept at £2 = 0.24. As expected, 
the two AK = 2 pairs feel a strong Y32 interaction so that the 
force causes large energy splittings in the diagram. As a con- 
sequence, an energy gap at proton number 98, and a one at 
106, show up. This result implies that nuclei in the mass- 
250 region, with proton number from 96 to 108, will possibly 
exhibit a strong triaxial-octupole effect with the Y32-shape. 
Particularly, isotopes of the element 98 (Californium) and 106 
(Seaborgium) can feel very strong nonaxial-octupole corre- 
lation. However, much less Y32 effect is observed around 
N=150 in the neutron single-particle diagram. 

It will be more convincing if the above arguments at the 
mean-field level are supported by some quantitative calcula- 
tions. If the observed 2~ bands are indeed a manifestation of 
the triaxial-octupole correlation, the calculations that properly 
include the octupole degree of freedom should be able to re- 
produce them. To demonstrate this, we calculate the low-lying 
spectrum in these nuclei by using the RASM [33]. TheRASM 
extends the Projected Shell Model l34tl applicability by break- 
ing the reflection symmetry in the basis and carrying out par- 
ity projection in addition to angular momentum projection. 
In this way, the octupole correlation is efficiently described 
and the model is thus particularly suitable to be applied to 
nuclei that show octupole effects. The trail wave function of 
the triaxial RASM is written in terms of parity- and angular- 
momentum-projected multi-qp states 



Kk 



(i) 



where the index K labels deformed basis states, and P p 
and PIj K are the parity-projection operator and angular- 
momentum-projection operator, respectively. For even-even 
nuclei, \<f> K ) in Eq. (Q]) is 
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where a v and a\ are the qp creation operator for neutrons and 
protons, respectively, with respect to the qp vacuum |0). The 
qp states are obtained from a deformed Nilsson calculation 
followed by a BCS calculation, in a model space with three 
major shells for each kind of nucleon (major shells 5, 6, 7 for 
neutrons and 4, 5, 6 for protons). The deformed Nilsson calcu- 
lation is performed with the quadrupole term and the octupole 
term; The former term breaks the rotational symmetry, and the 
later the reflection symmetry, in the qp states. The utility of 
symmetry breaking is to introduce efficiently the correspond- 
ing correlations in the basis states. The symmetries thus spon- 
taneously broken are recovered in the total wave function by 
projection (see Eq. (fl}). 

The two-body Hamiltonian has the form 

X=2 n=-X fi=-2 

(3) 
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In Eq. (O, Ho is the spherical single particle Hamilto- 
nian, the second term contains the quadmpole-quadrupole and 
octupole-octupole interactions, and the third and the last term 
are the monopole pairing and quadrupole pairing force, re- 
spectively. The interaction strengths in the second term are 
related to the deformation parameters through self-consistent 
relations 




5 h03 n (Q2o}„ + hO)p(Q2o)p' 
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7 tz(O ll (Q32}n+tl(O p (Q3 2 )p' 



(4) 
(5) 



Therefore, for given deformations, % are calculated, and thus 
not regarded as parameters. In the present work, we fix 
£2 = 0.235 for all nuclei; this value is verified to be consistent 
with that obtained from the TRS calculations. The octupole 
deformation parameter £32 influences sensitively the energy 
separation of the 2~ band from the ground state. We find that 
to reproduce data, the best suitable values are: £32 = 0.105 
for 246 Cm, 0.118 for 248 Cf, 0.110 for 250 Fm, and 0.107 for 
252 No. We note that 248 Cf with Z = 98 requires the largest £32 
among the isotones. The monopole pairing strength Go takes 
the usual form 



Go 



(6) 



where the minus (plus) sign stands for neutrons (protons), and 
g\ = 20.36 and g2 = 1 1 .26 are fixed for all the nuclei studied 
in this paper. The quadrupole pairing strength G2 is assumed 
to be proportional to Go, with the proportional constant 0. 13. 

Diagonalizing the Hamiltonian (f3]i in the space spanned by 
parity- and angular-momentum-projected multi-qp configura- 
tions (O, we obtain a set of energy levels. From them, we 
plot in Fig. 3 the ground band (the lowest A/ = 2 rotational 
sequence with positive-parity) and the 2~ band (the AI = 1 
sequence with negative-parity), both of which are compared 
with available data. As one can see in Fig. 3, agreement be- 
tween theory and experiment is excellent. We also obtain a 
good agreement with the unpublished 2~ band in 250 Fm 13511 . 
A careful look indicates that, to compare with the neighboring 
isotones, the nucleus 248 Cf has a notably low 2~ band. It is 
remarkable that this low 2~ band can be well reproduced by 
using an enhanced £32 in the calculation (see the employed £32 
values in the last paragraph), reflecting the strongest triaxial- 
octupole correlation in this nucleus. The physics behind can 
be understood by looking at the right part of Fig. 2. It can be 
clearly seen that there is an energy gap at the proton number 
98. The emergence of a gap is a consequence of the maxi- 
mum triaxial-octupole interaction, and therefore, 248 Cf shows 
the lowest 2~ band. For the same reason, we expect very low 
energy 2~ band in the Sg nuclei (the element 106). 

These low-lying 2~ bands are collective in nature. The de- 
gree of octupole collectivity is generally determined by an en- 
hancement of reduced E3 transition probability B(E3) over 



the single-particle estimate. The reduced rate of electromag- 
natic transitions between the eigenstates | V// p m), induced by a 
spherical tensor operator M^, is 



B(MX;i^f) 



2^+1 



2/i + l 

and the reduced matrix elements become 

1 



(7) 
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The calculated reduced transition probabilities from ground 
state 0+ to the 3" state of the 2" band, B(£3;0 + -> 3"), 
are shown in Fig. 3 in Weisskopf (single-particle) units 
(W.u.) (see the values marked by 'x'). Again, the most en- 
hanced triaxial-octupole collectivity appears in 248 Cf, where 
the B(E3) value is the largest. 

As we are exploring a completely unknown island of super- 
heavy elements, it is very useful to understand the structure of 
the mass-250 nuclei, which are the gateway to the superheavy 
region. Most studies on the excited configurations to date have 
been focused on the multi-quasiparticle structure (see, for ex- 
ample, Refs. l36l [13l). Recently, alternative parity bands in 
the heaviest nuclei are predicted B38I1 . The excited configura- 
tions are sensitive to the detailed shell structure, and thus can 
bring more useful information. The degeneracy of the pro- 
ton i'i 3 /2 and fj/2 states found in the single-particle diagram 
may be accidental; however, it creates the right condition to 
form the strong triaxial-octupole correlation that explains the 
experimental data. Thus, the present proposal seems to set up 
a new, strong constraint on all mean field calculations that are 
employed to predict shell structure for the superheavy mass 
region. 

In conclusion, the nonaxial-octupole F32 correlation has 
been found to be important in the transfermium and super- 
heavy region for the elements from 96 to 108, with the 
strongest effect in the elements 98 and 106. We have 
demonstrated that if this correlation is considered in realistic 
RASM calculations, the experimentally observed low-energy 
2~ bands in the AT = 150 isotonic chain can be easily de- 
scribed. The present work has thus suggested that this exotic 
nuclear effect, which has been discussed theoretically in other 
mass regions but with no convincing experimental evidences, 
may show up in the heaviest nuclei. The unique role of the 
nonaxial-octupole degrees of freedom gives weight to the re- 
cent prediction of tetrahedral symmetry in atomic nuclei. It 
has also set up a constraint on theoretical single-particle states 
for the superheavy region. The reduced E3 transition proba- 
bilities for the Coulomb excitation from the + ground states 
to the 3~ states of the 2~ bands are predicted to be strong. 
The future measurements of E3 transition rates would be cru- 
cial for the final determination of the nonaxial octuple nature 
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FIG. 3: The RASM calculations for the ground band and 2 band in the N = 150 isotones. The theoretical results are compared with available 
data. Note the very low 2~ band in 248 Cf. The calculated B(E3; 0+ — > 3~) value in Weisskopf unit is marked by 'x'. 



of the 2~ bands and, consequently, their relation to the tetra- 
hedral symmetry. 
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